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ABSTRACT 

The inclusive density of charged secondaries produced in high energy heavy ion col- 
lisions are calculated in the framework of multiple scattering theory and quark-gluon 
string model. The obtained results are in agreement with the data only assuming rather 
large inelastic screening/percolation effects. We predict the total saturation of relative 
inclusive densities in dependence on the number of interacting nucleon. Predictions for 
LHC are also given. 
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1. Introduction 



The processes of high energy heavy ion collisions can be considered in the framework of 
multiple scattering theory Many references can be found in reviews |], |2|, 0, f| • Contrary 
to the case of hadron-nucleus interactions, in the case of heavy ions it is impossible to 
sum analytically]] the contributions of all diagrams of Glauber approach. 

However, the contributions of the most important diagrams can be accounted for 
analytically, that allows one to calculate the integral cross sections of different processes, 
distributions on the number of interacting nucleons, multiplicities of secondaries and 
their inclusive distributions in reasonable agreement with the data. 

It is well-known that in high energy hadron-nucleus collision there exists inelastic 
screening || |7J which is confirmed experimentally, especially for the case of hadron- 
deuteron interactions. The same inelastic screening should exist in high energy heavy 
ion collision. This effect is very small for integrated cross sections (because many of them 
are determined by geometry), but it is very important || for the calculations of secondary 
multiplicities and inclusive densities. Similar results can be obtained §|, [10| in the frame- 



work of string fusion [ II], or percolation [12 models where string fusion/percolation 
effects directly correspond jl3) to the pomeron interactions. 

In the present paper we will give the theoretical estimations for relative inclusive 
densities, integrated multiplicities and rapidity distributions using both Gribov's reggeon 
diagram technique |L4| and percolation approach as well as Quark-Gluon String Model 
(QGSM) [0,0. 

We will also compare the above estimates with results from the Dual String Model 
(DSM) || [17j and find them remarkably similar, although the DSM is based on a different 
approach: it takes as input quantities measured in pp and/or pp collisions, and screening 
is realised via fusion of strings exchanged in partonic sea collisions, which can lead to 



percolation effects 18 



We will show that there exists a qualitative difference between secondary production at 
CERN SpS and RHIC energies and this difference can be described only by accounting for 
rather large inelastic screening/percolation/string fusion effects. These effects decrease 
the inclusive density about two times, that is in agreement with previous estimations || . 

For the distributions calculated in the present paper in DSM it is string fusion that 
leads to saturation with the number of participant nucleons, percolation phase transition 
is not felt. 

2. Spectra of secondaries and multiplicities in heavy ion collisions 

We start with standard multiple scattering theory consideration without inelastic 
screening. At high energies the spectrum of a secondary particle h produced in the 
central rapidity region of nucleus A - nucleus B collision is proportional to the same 



1 It can be done with the help of Monte Carlo calculation of multidimensional integral [^). 
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Figure 1: Diagrams corresponding to the rigid target approximation for nucleus A - 
nucleus B interaction. Only interacting nucleons are shown as the open points. 



spectrum in NN collision [pL9| 



da(AB -> hX) _ da(NN -> hX) 



dy 



dy 



However at existing energies the right-hand side of this expression is significantly 
smaller due to energy conservation corrections. The values of these corrections depend 
on the classes of multiple scattering diagrams which are taken into account. This problem 
was analysed in |20| and the results obtained in the so-called rigid target approximation 

4 GeV per nucleon) experimental 
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J were in agreement |22j with the low energy 
data on heavy ion collisions. 

Rigid target approximation corresponds to the account for the contributions of dia- 
grams shown in Fig. 1. Every nucleon of a beam nucleus A can interact several times 
with a target nucleus B, however every nucleon of B nucleus can interact only one time. 
This asymmetrical picture seems to be reasonable for nucleus A fragmentation, because 
the main part of shadow corrections is accounted for. For fragmentation of nucleus B 
the result seems to be inconsistent. So the inclusive spectrum of secondary h from AB 
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collision can be written |2(J as 



1 da(AB^hX) . . l j . Ar 1 da(NB^hX) 



/ s h / x 1 da(NA -> /iX) 
+ 9(-y)R h B (-y)<N B > — d 1 (2) 



where y is the rapidity of secondary ft, in c.m.s., < > and < Nb > are the averaged 
numbers of participating nucleons for every incident nuclei and the functions R\ B (y) 
account for the energy conservation effects. These last quantities can be found from the 
calculated inclusive spectra in NA and NB collisions, and they differ from unity about 
20% at CERN SpS energy and about 10% at RHIC energy 

Now we have to calculate the inclusive cross sections of secondary particle h in NA 
and NB collisions. It can be done with the help of the quark-gluon string model (QGSM) 
15 , [I~6|j . This model is based on the Dual Topological Unitarization (DTU) and it de- 



scribes quite reasonably many features of high energy production processes including the 
inclusive spectra of different secondary hadrons, their multiplicities, KNO-distributions, 
etc., both in hadron-nucleon and hadron-nucleus collisions [15, 16, 23]. 



High energy interactions are considered in QGSM as proceeding via the exchange 
of one or several pomerons and all elastic and inelastic processes result from cutting 
through or between pomerons ^4j. The possibility of exchanging a different number of 
pomerons introduces absorptive corrections to the cross sections which are in agreement 
with the experimental data on production of hadrons consisting of light quarks. 

Each pomeron corresponds to a cylindrical diagram, Fig. 2a, and thus, when cutting 
a pomeron two showers of secondaries are produced (Fig. 2b). The inclusive spectra 
of secondaries are determined by the convolution of diquark, valence and sea quark 
distributions u(x, n) in the incident particles and the fragmentation functions G(z) of 
quarks and diquarks into secondary hadrons. Both the initial quark distributions and 
the fragmentation functions of quarks and diquarks are constructed using the reggeon 
counting rules |^5| . 



The diquark and quark distribution functions depend on the number n of cut pomerons 
in the considered diagram. In the following we use the formalism of QGSM. In the case 
of a nucleon target the inclusive spectrum of a secondary hadron h (its rapidity, y, or 
Feynman-x distribution) has the form |T5[: 



dn 1 da x E da ^ h 

dy a inei dy a inei dx F ^ 

where the functions 4>n(x) determine the contribution of diagrams with n cut pomerons 
and w n is the probability of this process. Here we neglect the contributions of diffraction 
dissociation processes which are comparatively small in most of the processes considered 
below. It can be accounted for separately |TH . 
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Figure 2: Cylindrical diagram which corresponds to the one-pomeron exchange contribu- 
tion to elastic pp scattering (a) and its cut which determines the contribution to inelastic 
pp cross section (b). 
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For pp collisions 

4tp{x) = n)/*(s_, n) + n)/£(x_, n) + 2(n - (x +J *)/?(*-, n) , (4) 

x± = 4m^ /s + x 2 ± x] , (5) 

where f q and / s correspond to the contributions of diquarks, valence and sea quarks 
respectively. They are determined by the convolution of the diquark and quark distri- 
butions with the fragmentation functions, e.g., 

fg(x+,n) = / u q (x 1 ,n)Gg(x + /x 1 )dx 1 . (6) 

The diquark and quark distributions as well as the fragmentation functions are deter- 
mined from Regge intercepts. Their expressions are given in Appendix. 

The calculation of the inclusive spectra on nuclear targets is similar to Eq. (3). The 
only difference is that in the case of interaction, say, with three target nucleons and with 
exchange by five pomerons, it is necessary to account for all possible permutations. 

The multiplicities of the produced secondaries can be obtained by integrating the 
inclusive spectra. 



3. Comparison with data 

Many examples of the QGSM description of experimental data can be found in |15|, 
|16[ |23| , [26| , p7f . Now we are more interested in particle yields in central (mid-rapidity) 
region. The data on secondary ir and production in pp collisions at ISR energies |28|] 
at 90° in c.m.s. are presented in Figs. 3a and 3b. We slightly change the pomeron theory 
parameters, namely its intercept was assumed to be «p(0) = 1 + A with A = 0.09. This 
values of A results in slightly better description of the ISR data. It is practically the same 
as Donnachie-Landshoff p9f value A = 0.08 and rather smaller than the value A = 0.17 
which was obtained in PD[. The description of secondary protons and antiprotons can 



be found in |2q| . One can see that these data are described by QGSM quite reasonably. 

Now we can consider the secondary production in heavy ion collisions. The results of 
calculations of inclusive density of charged secondaries in mid-rapidity region for CERN 
SpS energy are in agreement with the data, whereas for RHIC energies are more than 
two times larger than the RHIC experimental data, and it is impossible to change QGSM 
parameters to have significantly better description for both pp and heavy ion cases. 

This disagreement is rather expected. At comparatively low energies (not higher 
than several GeV per nucleon) the heavy ion collisions can be considered with the help 
of standard Glauber approximation, and the inclusive spectrum of any secondary h 
produced in the central (mid-rapidity) region is described by the diagram shown in Fig. 
4a, i.e. by the contribution of a single nucleon-nucleon blob. This diagram immediately 
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Figure 3: Secondary negative (a) and positive (b) pion and kaon yields at ISR energies 
p8| at 90° in c.m.s. together with the data for all charged secondaries and with QGSM 
predictions. 

gives Eq. (1) for heavy ion inclusive cross section. The contributions of all other diagrams 



cancel each other due to AGK cutting rules [[24 



At high energies the new diagrams appear which include the interactions of pomerons. 
The simplest example of such diagrams is so-called triple-pomeron diagram which de- 
scribe the diffractive production of a large mass M jet in pp interactions. At low energies 
the contribution of such diagrams to hadron-nucleus or nucleus-nucleus interaction is 
suppressed by the longitudinal part of nuclear form factor [Q 

F A (tmin) ~ exp (R 2 A t min /3) (7) 

for Gaussian distribution of nuclear density, where 

tmin « (m N M 2 /s) 2 . (8) 

When energy becomes high enough, the value of t min becomes very small. So the dis- 
cussed contribution can be significant and namely these diagrams determine the inelastic 
screening corrections to hadron-nucleus cross sections. 

One example of a diagram with pomeron interaction for heavy ion interaction is shown 
in Fig. 4b. Contrary to hadron-nucleus case the contribution of such diagram can be 
estimated from the processes of high mass jet production in mid-rapidity region and 
with two large rapidity gaps, see for example Fig. 4c. The numerical contribution of all 
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Figure 4: Diagrams for inclusive cross sections for nucleus A-nucleus B collisions in 
Glauber approximation (a) and with accounting for the interactions of pomerons (shown 
by wave curves) (b). An example of inelastic processes of pp interactions which determine 
one of the needed vertices of the multipomeron interactions (c). 
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such diagrams is rather unclear because the number of diagrams is very large and the 
vertices of multipomeron interactions are unknown. For example, in || the Shwimmer 
model was used for the numerical estimations. However all such estimation are model 
dependent. 

The numerically significant contribution of the considered diagrams to inclusive den- 
sity is suppressed quadratically, by both longitudinal form factors, FA(t min ) and F B {t min ). 
So we can observe their influence at energies quadratically higher in comparison with 
the energies region where the inelastic screening effects are observed in hadron-nucleus 
scattering. One can see that the RHIC energies are of the needed order of magnitude. 

The reasonable possibility to estimate the contribution of the diagrams with pomeron 
interaction comes from percolation theory. We assume that if two or several pomerons are 
overlapping, they become one pomeron. It means that the inclusive density is saturated, 
it reaches its maximal value at given impact parameter. This approach has one free 
parameter - the critical number of pomerons in one squared fermi. Technically it is 
more simple to bound the maximal number of pomerons, n max , which can be emitted by 
one participating nucleon for the given pair of colliding nuclei. All model calculations 
become rather simple because above the critical value every additional pomeron cannot 
contribute to the inclusive spectrum. 

The Dual String Model on the other hand is based essentially on Dual Parton Model 



ideas |3l] with the inclusion of strings |TTJ which may interact in the transverse plane 
and fuse, thus reducing their contribution to the final state. Hadrons are considered as 
made up of constituent quarks (valence and sea quarks): the valence- valence diagram 
corresponds to single inelastic scattering and the wounded nucleon model [021, while the 
diagram (including gluons) corresponds to additional inelastic multiple scattering 
contributions. These contributions may be internal, parton multiple scattering within 
the original valence- valence contribution, or external, involving other nucleons. It should 
be noticed that the DSM applies only to symmetrical collisions (A=B). 

The main equation of the DSM was originally written for the rapidity particle density 

as 

JL = N A [2 + 2(k - l)a] h+(v- N A )2kah, (9) 

°y n a n a 

where h is the height of the rapidity plateau for valence- valence collisions, ah the height 
of the plateau, 2k the average number of strings produced in a nucleon-nucleon 

collision and N A is the number of participating nucleons from nucleus A; finally v is the 
average number of nucleon-nucleon collisions, which from elementary multiple scattering 
arguments satisfies at high energy v ~ N A 3 . Notice that the first term on the 
right-hand side is exactly A^i multiplied by the contribution of one nucleon-nucleon 
collision, while the second term results from type collisions. One should notice 

that the number of nucleon-nucleon collisions is N A + (y — Na) = v and the number 
of strings is N A [2 + 2(k — 1)] + [y — N A )2k = 2kv. Assuming that h and a are energy 
independent (constant plateaus) the energy dependence of dN/dy\ pp , taken in this model 
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Pb collisions at y/s = 17.3 
multiplied by 1/2) and for Au — Au at 
7|) and at y/s = 56 GeV per nucleon (open 



Figure 5: Relative inclusive densities of secondaries for Pb 
GeV per nucleon (black points, [[34], 
y/s = 130 GeV per nucleon (open points, 
triangle, f5B|). In Fig. 5a the solid curve shows the results of calculations accounting 
for the interactions of pomerons with n max = 1.67. Dashed curves present the results of 
the Glauber approximation (without percolation effects) for y/s = 17.3 GeV per nucleon 
(lower curve) and y/s = 130 GeV per nucleon (upper curve). Dotted curve shows our 
predictions for Au — Au at 56 GeV per nucleon with n max = 1.67, and dash-dotted curve 
1 and 2 for Pb — Pb collisions at 5.5 TeV per nucleon with n max = 1.67 and n max = 1, 
respectively. In Fig. 5b the results of DSM calculations at y/s = 130 GeV per nucleon 
with and without string fusion are shown by solid and dashed curves. 
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by interpolating experimental data, fixes the energy dependence of k, as explained in 
ref. |J, where it was found that a = 0.05 and h = 0.75; all results of the DSM shown in 
this paper are obtained with such values. 

By integrating Eq. (|9|) over the whole rapidity interval predictions are obtained from 
the DSM for the average multiplicity: the model's result is clearly linear but it overesti- 
mates the data. As argued in |T7| , the assumption that the plateau is proportional 
to the valence- valence plateau is not valid in the fragmentation region, effectively one can 
say that a = there. Fitting the experimental data allows us to estimate that Eq. (|9p 
is valid over the central region corresponding to 70% of total phase space. 

The results of calculations of the inclusive densities of the produced charged secon- 
daries per one pair on interacting nucleons dn c h/dy/(0.5N part ), where N part = Na with 
the help of Eq. (2) and QGSM are shown in Fig. 5a. We present the experimental data 
for Pb — Pb collisions at yfs = 17.3 GeV per nucleon |34|, [36j for \y\ < 1 and for 



Au — Au at yfs = 130 GeV per nucleon |!6], |37|] and at yfs = 56 GeV per nucleon for 



\rj\ < 10. One can see that CERN SpS data are described reasonably without any addi- 
tional screening (percolation). It was shown also in [Q and it seems to be natural because 
the suppression effects of t min discussed above should be rather large at comparatively 
small energy of CERN SpS. 

However, the same calculation at RHIC energy yfs = 130 GeV per nucleon gives the 
relative inclusive density two times larger (upper dashed curve in Fig. 5a) than the data, 
that is in numerical agreement with ||. The agreement with the data can be obtained 
only with suppression of multipomeron contributions, namely by using n max = 1.67 for 
every interacting nucleons. The results of DSM calculations with and without string 
fusion are shown in Fig. 5b, and they are in agreement with the curves in Fig. 5a. 

The same value of n max = 1.67 allows us to describe the PHOBOS point at 56 GeV 
(dotted curve in Fig. 5a). In the case of LHC energy yfs = 5.5 TeV per nucleon for 
Pb — Pb collisions with n max = 1.67 percolation effect decreases the relative inclusive 
density about 3 times (again similar to || prediction). The recent result of PHOBOS 
Coll. [[38] for central Au — Au collisions at yfs = 200 GeV per nucleon gives dn c h/dr] = 
650 ± 35 for \rj\ < 1 that is in agreement with our result dn c h/dr] « 600. 

However, the percolation effect at LHC energy can be even larger. The transverse 
range of a pomeron increases with incident energy as 



a' In s/ sq , (10) 



that is well-known experimentally from the shrinkage of slope parameter in elastic 
hadron-nucleon scattering. So in the percolation picture the value of n max at higher 
energies should decrease. The predictions for the relative inclusive density at LHC en- 
ergy with n max = 1 (maximal percolation) are shown in Fig. 5a by dash-dotted curve 2. 
Now the percolation effect decreases the relative inclusive density about 5 times. 

2 We accounted for the difference between rapidity and pseudorapidity distributions. 
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Let us note that Eq. (2) predicts independence of dn ch /dy/0.5NA on the Na value 
for Na larger than < Na >, that is confirmed by the dataQ. Really some rather weak 
dependence can exist because in central collisions every participating nucleon interacts 
more "centrally" on the average in comparison with minimum bias events. However 
this numerically small effect (~10-20% |3"9fl ) should be practically suppressed by perco- 
lation effects except of rather peripheral collisions, where the value of N part is small and 
percolation effects are not so important. 

The same calculations in the framework of DSM gives similar results as it is shown in 
Fig. 5b for PHENIX data at y/s = 130 GeV. Notice that contrary to the case of QGSM 
approach with Eq. (2) there is no saturation with the number of participants without 
string fusion. 

In Fig. 6 we present the calculated values of the relative inclusive density for secondary 
production in pp and Pb — Pb collisions as the functions of initial energy. One can see 
that due to percolation of pomerons emitted by different nucleons with similar impact 
parameters, the relative inclusive density in Pb — Pb collisions can be smaller than in 
the pp case. In the case of DSM model the growth with energy is more pronounced than 
in the QGSM due to the extrapolation of pp data with the second power of log s. 

Now let us consider the data ^Df on the total number of charged particles detected 
within range —5.4 < rj < 5.4, i.e. in the region which exclude only high fragmentation 
and diffraction regions. 

These data are in good agreement (Fig. 7) with our calculations in both QGSM and 
DSM including pp point. Our curves are very close to a straight line, again due to Eq. 
(2) for QGSM case. 

In |4(| it was observed also the difference in redistributions for the central and pe- 
ripheral Au — Au collisions. Qualitatively the same difference for rapidity distributions 
of charged secondaries produced in pp and Au — Au collisions at y^i =130 GeV per 
nucleon is presented in Fig. 8. Some shoulder in the dashed curve is connected with 
secondary proton contribution. It is an artifact of our approach, because we assume that 
all secondaries have the transverse momenta equal to their averaged values. 

The observed differences for rapidity distributions and their numerical values are 
rather important. They should be connected [41] with the possibility of multiple inter- 
actions of every nucleon in heavy ion collision. 



4. Conclusions 

Immediately after the first data of RHIC and their comparison with CERN SpS data 
we see for the first time in high energy physics the pomeron (secondary particle) density 
saturation (percolation effects or contribution from pomeron interactions). 

3 In the presented approach we can calculate inclusive density for Na smaller than < Na > only in 
DSM. 
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Figure 6: Predictions for the relative inclusive charged particles density for secondaries 
within range \y\ < 1 in pp (dashed curve, QGSM and dotted curve, DSM) and Pb — Pb 
(solid curve, QGSM and dash-dotted curves, DSM) collisions as the functions of initial 
energy. Two variants of solid curve correspond to percolation effects with n max = 1-67 
independently on energy and to its decrease from n max = 1.67 to n max = 1. 
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Figure 7: Total number of charged particles detected within range —5.4 < rj < 5.4 in 



Au — Au collisions at yfs = 130 GeV per nucleon [ 40| as a function of N part and its 
description by the QGSM with percolation effects and n max = 1.67 (solid curve) and 
DSM with string fusion (dashed curve). 
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Figure 8: Predicted by QGSM relative rapidity distributions of charged secondaries 
produced in Au — Au (solid curve) and pp (dashed curve) collisions at y/s =130 GeV 
per nucleon. 
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With accounting for this new effects with the help of simplest estimations and having 
only one free parameter (n max - the maximal number of pomerons emitted by one nu- 
cleon) we can reproduce many features, observed experimentally in heavy ion collisions. 
We predict the independence of relative inclusive density dn c h/dy / (0.5N part ) in the in- 
terval between minimum bias and central collisions, i.e. for A/2 < N part < 2 A p9|| , as it 
is shown in Fig. 5 for mid-rapidity events. We describe the same behaviour {n ch ~ N part , 
Fig. 6) for averaged multiplicities practically in all rapidity region. We reproduce the 
qualitative difference between the central and peripheral (or pp) heavy ion collisions. 

However the discussed new data show some problems for our approach. It seems 
the most important are rather small ratios of antibaryon to baryon in mid-rapidity 
region. For example, the ratio of p/p yield in the central region is about 0.6 whereas our 
approach predicts about 0.85 in agreement with calculations of [Kj. Possibly, the final 
state interactions |42| can change the situation. In the opposite case it means that we 
overestimate the contribution of sea quarks in the central region. 

We are grateful to N.Armesto for discussions. This paper was supported by grant 
NATO PSTCLG 977275. 



APPENDIX. Quark and diquark distributions and their fragmentation func- 
tions in QGSM 

In the present calculations we use quark and diquark distributions in the proton given 
by the correspondent Regge behaviour |15|, |25f . In the case of n-pomeron exchange the 
distributions of valence quarks and diquarks are softened than in the case of one-pomeron 
exchange due to the appearance of sea quark contribution. Also it is necessary to account 
that (i-quark distribution is more soft in comparison with u-quark one. There is some 
freedom |]16[ how to account for these effects. We use the simplest way and write diquark 



and quark distributions (for a>R = 0.5 and o>b = —0.5) as 

u uu (x,n) = C uu x a *- 2aB+1 (l-x)^ n -V- a * , (11) 
u ud {x,n) = C vd x a «- 2aB {l-x) n - 1 - a * , (12) 

(13) 

u u (x,n) = C u x~ aR (l-x) n+aR - 2aB+1 } (14) 

(15) 

u d (x,n) = C d x- aR (l-x)l {n ~ 1)+aR - 2aB+1 , (16) 

(17) 

Uu(x, n) = u 2 (x, n) = CuX~ aR 

X [(1 - x )n+a R -2a B -l _ & ^ _ ^n+2a B -2a B -l] , n > 1 , (18) 

u s (x,n) = C s x- aR {l-x) n+2aR - 2aB - 1 , n>\ , (19) 
where 5 = 0.4 is the relative probability to find a strange quark in the sea. The factors 
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Cj are determined from the normalization condition 
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Ui(x, n)dx = 1 (20) 



with sum rule 



/ y2 u i( x i ri)xdx = 1 . (21) 
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The fragmentation functions of quarks and diquarks into pions and kaons were changed 
a little in comparison with Refs. |15], |23[ to obtain a better agreement with the existing 
experimental data. We use the quark fragmentation functions in the form 

Gf = Gf = Gf = Gf = oo(l - zf- a fl - 0.7z) , (22) 

Gf = (l-z)Gf , (23) 

--a K (l-z) x , (24) 

G K U - = Gf = Gf = Gf = Gf = Gf = Gf (1 - z)(l - 0.95z) , (25) 

Gf = (l-z) Aa Gf , (26) 

z) x - aR (l -0.7z) , (27) 

1 - z)Gf , (28) 
with 

Act = a p - a = 1/2 , A = 2a' < p 2 >= 0.5 , a = 0.73 , a K = 0.24 . (29) 
Diquark fragmentation functions have the form : 

Gf = a (l - z) x+aR - 2a fl + 2z) , Gf = a Q (l - z )A+«*-2«b+1(i _ QQ ^ (3Q) 

Gf = Gf d = a (l - z) x+a *- 2a fl - 0.9*) , (31) 

Gf + =a K (l- z)W°*-°°) , Gf = Gf 1 + (1 2 ~ Z? , (32) 

Gf = Gf = a K (l - z f+z(«n-<*B) + i {l _ . 95z) . (33) 

The probability for a process to have n cut pomerons was calculated using the quasi- 
eikonal approximation |15], 

°° Op n_1 z k 

w n = (Jul <? n , a n = — (1 - e~ z tt) , (34) 

n=l nZ k=0 kl 

z = , °p = 8^7e A? , e = Hs/l GeV 2 ) , (35) 



with parameters 



A = 0.09 , a' = 0.21 GeV" 2 , 7 PP = 2.2 GeV" 2 



R 2 pp = 3.18 GeV" 2 , C pp = 1.35 
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